Hole spin relaxation in semiconductor quantum dots 
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Hole spin relaxation time due to the hole-acoustic phonon scattering in GaAs quantum dots 
confined in quantum wells along (001) and (111) directions is studied after the exact diagonalization 
of Luttinger Hamiltonian. Different effects such as strain, magnetic field, quantum dot diameter, 
quantum well width and the temperature on the spin relaxation time are investigated thoroughly. 
Many features which are quite different from the electron spin relaxation in quantum dots and 
quantum wells are presented with the underlying physics elaborated. 
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I. INTRODUCTION 



Recently, considerable interests have been devoted to 
spin-related phenomena in semiconductors due to the 
enormous potential of the spintronic devices. ^'^ Among 
these, properties of electron spins confined in semicon- 
ductor quantum dots (QD's) are essential to the pro- 
posed qubits in quantum computers and have there- 
fore caused much attention. '^■'''^'^'^'^ Many works calcu- 
lated the spin relaxation time (SRT) of electrons due to 
the spin-orbit coupling induced spin-flip electron-phonon 
scattering at very low temperatures, ^'^'^'^ where the dom- 
inant electron-phonon scattering arises from the piezo- 
electric potential. These works are based on perturba- 
tion theory where the spin-orbit coupling is treated as a 
perturbation in the Hilbert space spanned by Hq which 
does not include the spin-orbit coupling. Moreover only 
the lowest few energy levels of Hq are included in the 
theory. Recently we have shown that the perturbation 
method is inadequate in accounting for the electron struc- 
ture and therefore the SRT in semiconductor QD's: The 
SRT obtained from the perturbation approach used in 
the literature^'^'*'^ is several orders of magnitude smaller 
than the exact value. Very recently Wood et al. in- 
vestigated the SRT of a hole in QD's.^^ Again perturba- 
tion method is used and only the SRT induced by the 
electron-phonon scattering due to the deformation po- 
tential is considered. Many other effects such as strain, 
quantum well orientation and multi-subband effects as 
well as the effect from the electron-phonon scattering due 
to the piezoelectric coupling, have not been studied in 
their work. 

In the present paper, we investigate the hole SRT of 
GaAs QD's confined in quantum wells along (001) and 
(111) directions by parabolic potentials with strain in- 
cluded by exactly diagonalizing the hole Hamiltonian. 
We calculate the hole SRT due to the scattering with 
acoustic phonons by the Fermi golden rule after getting 
the hole energy spectra and the wavefunctions from the 
exact diagonalization. We discuss how the strain, QD 
radius, magnetic field, temperature and quantum well 



width affect the SRT. We show that strain on quantum 
wells of different growth directions affects the QD spin 
relaxation in totally different ways: for QD's in (001) 
quantum well strain changes the relative position of en- 
ergy levels of heavy hole and light hole; but for those in 
(111) quantum well, strain makes additional spin mix- 
ing and induces additional spin relaxation. Also we show 
that unlike the case of electrons where the SRT is mainly 
determined by the electron-phonon scattering due to the 
piezoelectric interactions, for holes both the hole-phonon 
coupling due to piezoelectric interaction and that due to 
deformation potential make important contributions to 
the spin relaxation process, although their relative im- 
portance changes under different conditions. 

The paper is organized as follows: In Sec. II we set 
up our model and Hamiltonian. In Sec. Ill we present 
our numerical results. We discuss the SRT of QD's in 
(001) quantum well in Sec. III(A). We first discuss a sim- 
ple case: a small QD without strain where we compare 
our results with those obtained from the perturbation 
method. Then we discuss strain dependence of the SRT 
when the confinement of the quantum well is very strong 
and there is only one subband. We finish Sec. III(A) by 
showing the strain, magnetic field, and QD radius depen- 
dence of the SRT in the case of large well width (multi- 
subband effects). Then we turn to the case of QD's in 
(111) quantum well in Sec. III(B). In Sec. III(C) we show 
the well width dependence of the SRT in both (001) and 
(111) quantum wells. We conclude in Sec. IV. 



II. MODEL AND HAMILTONIAN 

We use a simplified model to study the spin relaxation 
in QD's which are confined by parabolic potentials Vc(r) 
in a quantum well of width a. Due to the confinement of 
the quantum well, the momentum states along z axis are 
quantized. With the hard-wall approximation, the hole 
momentum states along z-axis are therefore characterized 
by the subband index . The total Hamiltonian is given 



2 



by 

H = Hh + H strain + Hph + Hint , (1) 

in which Hh is the 4x4 Luttinger Hamiltonian for holes. 
When the growth direction of the quantum well is along 



the (001) direction (z-axis) and the matrix elements are 
arranged in the order of = +|, — ^ and — |, Hfi 
can be written as^^ 
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In these equations, TOq denotes free electron mass; 71, 72, 
73 and K are Luttinger coefhcients; and and n'^ rep- 
resent the subband indices, w^"^ and cij°°^ in the two di- 
mensional confinement potential 14 (r) [Eq. (3)] represent 
the confinements experienced by the heavy hole and light 
hole respectively and arc given by u;^^'^ = h/ (m^u^cP) 



and 



,001 



m 



001 
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with m'^^^ 
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m^/ili + 72) and 



'mo/ ill — 72) standing for the effective masses 
of heavy hole and light hole in the direction perpen- 
dicular to the growth (001) direction and d represent- 
ing the QD diameter. By applying a magnetic field B 
along the growth (z) direction of the quantum well and 
adopting the Coulomb gauge A = (— 0), one has 



Px = {hk^ + ^) and Py 

From the Luttinger Hamiltonian Eq. (2) one can see 

that when the well width a is sufficiently small and only 
the lowest subband in QD is important, 5 = and -|-| 
(— |) states can only mix with — i (-1-^) states. Therefore 
there is no mixing between the spin- up and -down (±|) 
heavy hole states and between the spin-up and -down 
(±i) light hole states. Spin mixing between the spin- up 
heavy hole and the spin-down light hole or vice versa is 
negligible as the energy difference between the heavy hole 



and light hole states is usually too large. Nevertheless for 
larger well width where higher subbands are needed, S no 
longer equals to zero and the spin-up and -down heavy- 
hole states arc mixed with each other, mediated by the 
light-hole states. The same is true for the spin-up and 
-down light-hole states. 

The strain Hamiltonian H strain given by the Bir-Pikus 
Hamiltonian^'* has the form: 
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with the matrix elements being 

F = -{D^ + ^)Tr(.) + ^^e.. , 

G = -(D„-:^)TV(e)-^e.., 
H = Dd{ezx - iezy) , 



(7) 



xy 



(8) 

(9) 
(10) 

(11) 



Here Da, and Dd are the deformation potential con- 
stants, e is the strain tensor with eij denoting the tensor 
components. For (OOl)-orientcd zinc blende crystal the 
strain tensor components are given by*"'' 



^001 



^001 
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(12) 

(13) 
(14) 
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where ai and 02 are the lattice constants of epilayer 
(GaAs) and substrate materials, and C'u and C12 are the 
stiffness constants. One can see that for (OOl)-oriented 
zinc blende crystal, H = I = and the strain Hamilto- 
nian Eq. (7) has only the diagonal terms. Therefore the 
strain does not induce any extra spin mixing, but adjusts 



the relative positions of the heavy-hole and light-hole en- 
ergy levels. 

When the growth direction of the quantum well (z- 
axis) is along (111) direction, the hole Luttinger Hamil- 
tonian Hh is the same as that in Eq. (2), but with the 
matrix elements being replaced by^^ 



J 



P ± Q = (71 ± [P, + Py] + T 273)^ 

S = f (72 - 73)[P. + iPy? - ^(272 +73)- )2 _ (,^).] 



(1 - Sn;„n'J[Px - iPy] ' 



(15) 
(16) 

(17) 
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in Eq. (3) should be 



replaced by 

,111 _ fc //'^lll j2 



by cji" = h/{m.lYd'^) and w/" = ^/(m 



and mj^jp which are given 



^2) 111 _ 

;|| " ^ ^h\\ — 
77x0/(71+73) andmj^ii""'" = mo/(7i— 7,3) respectively. From 
Eqs. (15-17), one finds that differing from the previous 
(001) case, here S and R are nonzero even for the single 
subband case. This means when the growth direction 
is along (111) crystal direction, there is always mixing 
between the spin-up and -down heavy-hole states and 
the spin-up and -down light-hole states. 

The strain Hamiltonian Hstrain for (lll)-oriented zinc 
blende crystal is same as that in Eqs. (7-11), but now the 
strain tensor components arc given by-'^^ 



xy 



,111 



yy 



All 
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a = 



^yz ^zx 2 

Cii + 3Ci2 + 4C44 
2Cii + 4^12 - 4C44 



^(2-l/^"^)e|| 



(18) 
(19) 

(20) 



with C44 the stiffness constant. Substituting Eqs. (18- 
20) into Eqs. (7-11), one gets F = G and the off-diagonal 
elements H and / are no longer equal to zero. Conse- 
quently when one applies a strain on (lll)-oriented zinc 
blende quantum well, the change of energy levels of the 
light and heavy holes are almost the same and the strain 
introduces additional spin mixing between the spin-up 
and -down heavy/light holes. 

Hph in Eq. (1) is the Hamiltonian of acoustic phonons 

and is given by Hp^ = J2qx^Q^(^qx(^Q>^ with ^QA 
standing for the phonon energy spectrum of branch A 
and momentum Q. Two different hole- phonon scatter- 
ing mechanisms contribute to the spin relaxation for the 
temperatures we consider here. One is hole-phonon scat- 
tering due to piezoelectric coupling which is given by 

HZ = E ^Qa(«^-qa + «Qa) exp(iQ • r) (21) 

QA 



with Mqx being the scattering the matrix ele- 
ments. For longitudinal acoustic phonons Mq^i = 

327r e e^^^^ (3Qa:QyQ^) £ ^^^^ transvcrse acoustic 

phonons E,=i,2^^,*, = W^lQlQl + QlQl + 
^2q2 _ (3Q.Q»Q.)" j^ jjgj.^ ^ ^j^^ Qj^^g volume den- 



sity, ei4 represents the piezoelectric coupling constant 
and e denotes the static dielectric constant. The acous- 
tic phonon energy spectra wqa are given by ujqi = VsiQ 
for the longitudinal mode and ujQpt = VstQ for the trans- 
verse mode with Vgi and Vst standing for the correspond- 
ing sound velocities. Q = ^Ql + Ql + Ql- It is noted 

that this kind of scattering does not flip hole spin and 
therefore only when the hole wave function itself contains 
spin mixing can H^t contributes to spin relaxation. The 
other is hole-phonon scattering due to the deformation 
potential Hf^l . Hf^l can be derived from H strain Eq. 
(7) by substituting e-ij in Eqs. (8-11) by e-ij = (-ij+^ij, and 
splitting Hstrain into two parts: one contains all terms 
proportional to e?- and the other contains all terms pro- 
portional to e'ij-^^ The second part is therefore Hf^^ if 
e?- represent the strain tensor components caused by the 
sample and e'^j are the tensor components caused by the 
latticx; vibrations, e'^j can further be written in terms of 
normal-mode coordinates as 



QA 



2 y 2pVuJcix 



xiVixQj+VjxQiy^" 



(22) 



with p being the mass density of the material, Q stand- 
ing for the phonon wave vector and fjix representing 
the unit vector of polarization of A-phonon along i- 
dircction. For longitudinal mode rji =- Qi/Q and for two 



transverse modes rfti 
rjh = iQy,-Qx,0)/Q± with Q± = 



Qz,-QI)/QQl and 
Q2+g2.ii It is 
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pointed out that Hf^J is spin-flip scattering as it con- 
tains non-zero off-diagonal parts H and /, and therefore 
can cause spin relaxation even for the case when there 
is no spin mixing in the wavefunctions. Accounting for 
these two hole-phonon scattering mechanisms, one has 



TTpie 



H 



def 



int ^^int ^ ^^int 

We diagonalize the Luttinger Hamiltonian Hh in the 
Hilbert space \n,l,nz,cr) constructed by Hq which is 
taken to be the diagonal part of Hh: 

1*^)= E Cii^Jn,l,n,,<j) . (23) 

nlrizO' 



lA^QAp = |MQA(/|e^'^-'"|i)|2. It is noted that only 
when the eigenstates \i) and |/) contain spin mixing can 
■A^QA 7^ be possible. For hole-phonon scattering due to 
the deformation potential, |A1qaP = \{f\HQ\\i)\^ with 
Hq\ the matrix for the hole-deformation potential. As 
iJqA itself contains spin mixing, therefore it is not neces- 
sary to have spin-mixed initial and final states to ensure 
■A4qa 7^ 0. The total SRT r can be written as 



7 = E/^Er.- 



(28) 



Here HQ\n,l,nz,(T) = E^j Jn,l,nz,a) with 



{r\n,l,n^,a} 



2 . ,nzir 



X\l - sin( z) 
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mo 



(24) 



m 



[hn{2n+\l\ + l)-hujBl] 



h\\ 

, iheBn mo h^TT^ni 
± — \ —- ^ , (25) 



2mo TO^^ 2moa2 



^[hn{2n+\l\ + l)-hujBl] 
^1 

, heBK mn fi^TT^nl 

± 1 — . 

2mo rnf^ 2moa'^ 



(26) 



In these equations n = 0, 1, 2, • • • and I = 0, ±1, ±2, • • • 
are quantum numbers; ^ denotes the growth direction 



and 



Hz 



stand for the 



which can be (001) or (111); 
effective masses of heavy and light holes in the z-direction 
which are given by mo/m^°^ = 7i — 272, mo/mj^^ = 
71 - 273, m o/m^°^ = 71 -I- 272 and mo/m\l^ = 71 + 
273; O = -^/wq + UJ% with ujQ = h/{mo(P) and lob = 

eB/(2mo); Nr,,i = { ^^^+ 1 1 1 )') ^ "^^^^ " = moQ/h. 4' is 
the generalized Laguerre polynomial. 

The eigenfunction | ^E'^) in Eq. (23) is a mixture of four 
different components: spin-up and -down heavy-hole and 
light-hole states. We assign an eigenstate i to be spin up 
if the spin-up components are larger than the spin-down 
ones. An hole at initial state i with energy and a spin 
polarization can be scattered by the phonon into another 
state / with energy e/ and the opposite spin polarization. 
The rate of such scattering can be described by the Fermi 
golden rule: 



^ QA 

+ (^QA + l)'5(e/ -ei+ cjqa)] , 



(27) 



with uqx representing the Bose distribution of phonon 
with mode A and momentum Q at temperature T and 

Mq\ being the corresponding matrix elements. For 
hole-phonon scattering due to the piezoelectric coupling. 



in which fi = Cexp[— ei/(fcsT)] denotes the Maxwell 
distribution of the i-th level with C being a constant. 



III. NUMERICAL RESULTS 

It is seen from our previous discussion that hole spin 
relaxation in QD's is very complicated and is affected 
by many effects. When the quantum well is along (001) 
direction, for small well width (with only the lowest sub- 
band) only hole-phonon scattering due to deformation 
potential contributes to the spin relaxation; for large well 
width (with multi subbands) hole-phonon scattering due 
to both deformation potential and piezoelectric coupling 
contributes to the spin relaxation. Strain itself in this 
case cannot bring additional spin relaxation but influ- 
ences it by changing or even reversing the relative posi- 
tions of energy levels of heavy hole and light hole. Es- 
pecially when a minus strain makes the energy levels of 
heavy hole and light hole very close to each other, spin 
mixing between spin-up (-down) heavy hole and spin- 
down (-up) light hole cannot be neglected anymore and 
therefore the hole-phonon scattering due to piezoelectric 
coupling may contribute to spin relaxation also. Never- 
theless, when the quantum well is along (111) direction, 
regardless of the well width, there exists spin mixing be- 
tween spin-up and -down heavy (light) holes in the eigen 
functions of the Luttinger Hamiltonian. Therefore hole- 
phonon scattering due to both deformation potential and 
piezoelectric coupling contributes to the spin relaxation. 
Moreover, strain itself in this case makes additional spin 
mixing and induces additional spin relaxation. In this 
section we perform a comprehensive investigation to find 
out the relative importance of above mentioned effects 
under various conditions such as temperature, magnetic 
field, QD radius and quantum well width. The parame- 
ters used in the computation are listed in Table 1.1^.16,17 

In order to ensure the convergence of the energy spec- 
tra €£, we use sufficient basis functions to diagonalize the 
hole Hamiltonian Hh- For example, in a QD with B = 1 
T, a = 5 nm, d = 20 nm and without strain one has to 
use 100 basis functions to converge the lowest 40 levels; 
Nevertheless, when a = 20 nm one has to use 484 basis 
functions to converge the same levels. 
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P 

Vst 
Vsl 

ei4 
Cii 



5.3 X 10^ kg/m^ 
2.48 X 10^ m/s 
5.29 X 10^ m/s 

I. 41 X 10^ V/m 

II. 81 



e 

Da 

Dd 

Cl2 



12.9 

-6.7 eV 
-1.7 eV 
-4.55 eV 
5.32 



71 

72 
73 
K 

C44 



6.85 

2.1 

2.9 

1.2 

5.94 



TABLE I: Parameters used in the calculation. 
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FIG. 1: SRT vs. QD diameter d. Curve with •: exact 
diagonalization result with the energy sufficiently converged; 
Curve with X: Perturbation result; Curve with Q ■ exact 
diagonalization result but with only the lowest two heavy hole 
and the lowest two light hole levels used as basis functions. 
Curve with ■: exact diagonalization result but with only the 
16 energy levels of i?o given in the text as basis functions, 
(a): T = 0.2 K and (b): T = 4 K . 



A. QD's in (001) quantum well 

1. Small well width without strain 

We first consider QD's in a small (001) quantum well 
(a = 5 nm) without any strain where the lowest eigen 
states of Hq arc heavy holes and the separation between 
the heave and light holes is around 0.1 eV. Due to the 
small well width, only the lowest subband is needed in 
the calculation. Therefore, as pointed out above, only the 
hole-phonon scattering due to the deformation potential 
contributes to the spin relaxation. 

In Fig. 1 we plot the SRT as a function of the QD di- 
ameter for two temperatures. _B = 1 T in the calculation. 
Curves with • are the results obtained by the exact di- 
agonalization method with the energy levels sufficiently 
converged. Unexpectedly, differing from intuition as well 
as the results of an electron spin in QD's,^" it shows in 
Fig. 1(a) that the SRT increases with the QD diameter. 

To understand/check this result, we compare the re- 
sults from the exact diagonalization method with those 
from the perturbation approach widely used in the 
literature, ^''''^'^ but with the proper modification by in- 
cluding the second order corrections to the energy spec- 
trum as pointed out in our previous work.^" Here we 



treat the off-diagonal part of Luttinger Hamiltonian Hh 
as perturbation and calculate the SRT between the low- 
est two energy levels composed by two lowest heavy hole 
and two lowest light hole states of Hq: |0,0, l,cr) with 
cr = ±1 and The wave functions can be written as: 

(r|*T) = (r|0,0,l,^)-^(r|0,0,l,-^) , 



(r|f^) = (r|0,0,l,-^)-B(r|0,0,l,l) 



in which 



^= - 
B = - 



472?i2(7r2/a2 - a) - AheBn ' 
472/12(7^2 /a2 -a)+ AheBK 



(29) 

(30) 
(31) 



With the second order correction to the energy included, 
the energy difference between the spin- up I'^i) and spin- 
down l^fi) states can be written as: 



AE 



GheBn 
2mo 

GheBn 



+ l-4P(-£^o,o,i,| - -^0,0,1,1- 



E, 



0,0,1, 



3K7|(fcB)3/(4mo) 



2mo ['y2h'^{7r^/a^-a)]^-iheBKy 



(32) 



in which the first term represents Zeeman splitting. The 
SRT T can therefore be written as: 



Bf 



X exp(-|^)/2(g,) , 



(33) 



with Q = AE/{hvs), q ~ Qs'mO, — Qcos9 and 
HQz) = Stt^ sm{aqz/2)/{aqz['iTT'^ — (ag^)^]}. A stands 
for the branch of phonon: for longitude mode )Ci{6) = 
sin^ (6) cos^ (6) and for two transverse modes JCti{9) = 
isin(6')cos2(20) and /C(2(6') = i sin(6') cos2(6') respec- 
tively. The SRT's calculated from Eq. (33) are plotted 
as the curve with x in Fig. 1(a) which coincides with 
the curve with O obtained from the exact diagonaliza- 
tion method but with exactly the same four basis func- 
tions used in the perturbation method as basis. One can 
see that the SRT r does increase with the QD diame- 
ter at T = 0.2 K. From Eqs. (30-33) one finds that the 
SRT depends on the diameter only through a, which can 
be approximated into a = l/dP in the case wq ^ ojb- 
Moreover, the mixture of wave functions and the energy 

difference AE hardly change with d for tP" ja} ^ l/d^. 

2 

Therefore only the exponential term exp(— |^) in Eq. 
(33) decreases with d. As a result larger QD diameter 
corresponds to longer SRT at low temperature. 

As shown in our previous work that the right perturba- 
tion approach (with the second order corrections to en- 
ergy spectrum included) with the lowest few levels of Hq 
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as basis functions may lead to totally opposite trend from 
the exact diagonalization method with sufficient number 
of basis functions. In order to rule out this possibility in 
the present analysis, we plot in Fig. 1(a) the SRT calcu- 
lated from exact diagonalization method but with sixteen 
eigen functions of Hq as basis functions, i.e., |0,0, 
and |0, 1, 1, (t) with a = ±^ and ±|. It is seen from the 
figure that it produces the same T-d dependence. 

In Fig. 1(b) we plot the same curves but at T = 4 K. 
It is seen that at this temperature the SRT has a max- 
imum as a function of the diameter d. This is because 
at high temperatures the scattering between the higher 
energy levels becomes important. These high energy lev- 
els are arrayed very close to each other. The increase of 
diameters makes more levels into the scattering channels 
and thus induces a faster spin relaxation if it overcomes 
the opposite tendency described above. The perturba- 
tion approach only includes the lowest two heavy-hole 
and two light-hole energy levels and therefore can not 
get the maximum feature here. 
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FIG. 2: SRT vs. temperature T for QD with o = 5 nm 
and d = 20 nm. (a): Results under different magnetic fields. 
Curve with m: B = 0.2 T; Curve with k: B = 0.6 T; Curve 
with •: B = 1 T. (b): Results at B = 1 T with contributions 
from different branches of piionons specified. Curve with ■: 
Contribution from longitudinal phonons; Curve with Con- 
tribution from the transverse phonons of first branch; Curve 
with Contributions from the transverse phonons of the 
second branch; Curve with •: The total SRT. 

In Fig. 2(a) we plot the SRT as a function of the tem- 
perature for a QD with a = 5 nm and d = 20 nm under 
three different magnetic fields. From the figure one finds 
that the SRT decreases with the temperature as higher 
temperature leads to larger number of phonons ngx and 
consequently a larger transition probability. Moreover, 
as pointed out in our previous work,^° smaller magnetic 
field makes the SRT decrease faster with the tempera- 
ture. This is because the energy intervals between dif- 
ferent energy levels are small in the presence of a small 
magnetic field. This leads to a faster response to the 
temperature. In Fig. 2(b) we further specify the contri- 
butions from different branches of phonons. It is stressed 



again that only the holc-phonon scattering due to the 
deformation potential contributes to the spin relaxation 
here. As the temperature is below 4 K, the scattering 
from spin-up to spin-down states, in which phonons are 
emitted, is much larger than that from the reverse pro- 
cess, in which phonons are absorbed. Therefore, unless 
specified, the scattering rate i is defined to be the scat- 
tering from spin-up to spin-down states throughout this 
paper. It is seen from the figure that the SRT here is 
determined by the transverse modes. 

2. Small width with strain 
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FIG. 3: SRT versus strain for a QD in (001) quantum well at 
a = 5 nm, d = 20 nm and B = lTatT = 4K. Solid curve: 
total SRT; Chain curve: SRT induced by the holc-phonon 
scattering due to piezoelectric coupling; Dashed curve: SRT 
induced by the hole-phonon scattering due to deformation 
potential. 

As pointed out in the previous section, strain on (001) 

quantum well can change or even reverse the relative po- 
sitions of energy levels of heavy and light holes. Now we 
turn to investigate QD's under different strains in (GDI) 
quantum well with a = 5 nm, d = 20 nm and B = 1 
T at T = 4 K. The strain is adjusted by changing the 
strain tensor component e|| in Eqs. (12) and (13). We 
plot the SRT versus in Fig. 3 with eji obtained by 

substituting the lattice constants of GaAs and InAs for 
a2 and oi respectively in Eq. (12). It is seen from the fig- 
ure that when > 0, the SRT r increases with applied 
strain. This is because the positive strain enhances the 
gap between the heavy hole and light hole. Nevertheless, 
the flip from the spin-up heavy hole to the spin-down 
one is determined by the light-hole components in the 
wavefunctions. Increasing the gap between the heavy 
and light holes greatly reduces the spin relaxation and 
leads to the increase of SRT. When we apply a negative 
strain, as the gap decreases, the SRT decreases as shown 
in the figure. Particularly at e||/e|| = —0.3 the lowest 
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two energy states change from the heavy-hole states to 
the hght-hole ones and the SRT shows a minimum. When 
e||/e|j < —0.3 the initial spin states are light holes and in- 
creasing strain along the negative direction enhances the 
gap of light hole and heavy hole again and therefore the 
SRT is enhanced again. 

SRT around — —0.3 needs more address. Around 
this point, the energy levels of heavy hole and light hole 
are close to each other and therefore as said before that 
the hole-phonon scattering due to piezoelectric coupling 
is able to contribute to the spin relaxation. This can be 
seen in the inset of Fig. 3 where SRT due to the piezo- 
electric coupling is plotted as a chain curve and that due 
to deformation potential is plotted as dashed one. The 
sohd curve is the total SRT. Very close to e||/e|| = —0.3, 
the SRT is determined by the hole-phonon scattering due 
to piezoelectric coupling. However, the contribution of 
piezoelectric coupling decays dramatically with little de- 
viation of the strain from —0.3. When e||/ejj < —0.29 
or > -0.284, the SRT is totally determined by the hole- 
phonon scattering due to the deformation potential. 



3. Large width 




-1.0 -0.5 0.0 0.5 1.0 

^11 AS 

FIG. 4: SRT vs. strain for a QD in (001) quantum well at 
a = 20 nm, d = 20 nm and B = 1 T at T = 4 K. Solid curve: 
total SRT; Chain curve: SRT induced by the hole-phonon 
scattering due to piezoelectric coupling; Dashed curve: SRT 
induced by the hole-phonon scattering due to deformation 
potential. 

Now we turn to investigate the SRT in a QD of a = 20 
nm, d = 20 nm and B = 1 T at T = 4 K. For such a well 
width, one has to include states with > 2. Prom Eqs. 
(2-5) one can see that the scattering between different 
subbands makes ±| states mix with ±5 states. There- 
fore, hole-phonon scattering due to the piezoelectric cou- 
pling makes contribution to the spin-flip scattering with 
or without strain. 



In Fig. 4, we plot the SRT as hmction of eii/ejj. It 
is seen from the figure that the SRT shows a minimum 
around e||/e|| ~ —0.01. This is because for large quan- 
tum well, the lift of the F-point degeneracy is very small 
and the lowest heavy-hole and light-hole states are very 
close to each other. Again, the SRT increases with the 
applied positive/negative strain due to the separation of 
the heavy hole and light hole states. Differing from the 
case of single subband, here the hole-phonon scattering 
due to the piezoelectric coupling makes strong contribu- 
tion to the SRT even without strain. Nonetheless it is 
shown in the figure that when the strain is slightly 

deviated from —0.01 the SRT is mainly determined by 
the hole-phonon scattering due to deformation potential. 



T 1 1 r 




J I L 



30 35 40 45 50 
d (nm) 

FIG. 5: SRT vs. the diameter d for QD's at a = 20 nm and 
B = 1 T under different strains. Curve with •: ei|/e|| = 0; 
x: eii/ej} = 0.07; □: eii/ejj = -0.07; A: e\\/e\ = -0.09; 
4:e||/e^ = -0.3. (a) T = 0.2 K; (b) T = 2 K. 

In Fig. 5 we plot the SRT as a function of QD di- 
ameter d when a = 20 nm under different strains with 
solid curves for the strain-free case, chain curves for pos- 
itive strain cases and dashed curves for negative strain 
cases. Similar to the case of small well width without 
strain in Fig. 1, for large well width without strain here 
the SRT also increases with the diameter for low tem- 
perature (T = 0.2 K) [Fig. 5(a)] and shows a maximum 
for higher temperature (T = 2 K) [Fig. 5(b)]. Strains 
keep these trends. Nevertheless for small negative strain 
—0.2 < e||/e|j < the variations become smoother. This 
feature can be understood as follows: For negative strain 
the heavy hole states intercept with the light hole ones 
at high energies. The states around these intercepting 
points are particular efficient in spin relaxation. As said 
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in the previous section, the high energy states are ar- 
rayed very close to each other. The increase of diameter 
drives more states into the scattering channel and partly 
compensates the tendency of decrease. For larger neg- 
ative strain (e||/e|| < —0.2), the heavy and light holes 
arc separated again and the change of the SRT with the 
diameter becomes fast again as shown in the figure. 



10= 



10'^ 



10^ 
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+ _ 




0- o- o 
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B(T) 

FIG. 6: SRT vs. the magnetic field B for a QD in (001) 
quantum well at a = 20 nm, d = 20 T and T = 4 K under 
different strains. Curve with •: e||/e|| = 0; ■: e||/e|| = 0.01; 
♦ : eii/e;; = 0.03; Q: en/ejj = -0.01; □: e|i/eU = -0.02; x: 
e||/£;| = -0.03; +: en/ej} = -0.2. 

The SRT as a function of magnetic field B under dif- 
ferent strains at T = 4 K is plotted in Fig. 6. As shown in 
Eqs. (30) and (31), the spin mixing is enhanced with the 
increase of the magnetic field. Therefore r decreases with 
B. Moreover, similar to the case of the diameter depen- 
dence of the SRT, one finds that the t-B dependence 
becomes weak when a negative strain is in the range 
—0.2 < f:||/e|| < 0. This is because when the heavy-hole 
states intercept with the light-hole ones, Zeeman splitting 
which appears in the denominators of the spin mixing co- 
efficients [e.g., in the denominators of Eqs. (30) and (31)] 
becomes important which partially compensates the in- 
crease of the spin mixing with the magnetic field {ie., 
B in numerators of spin mixing coefficients). Therefore, 
the SRT changes with B slowly. However, for a larger 
negative strain which makes a big separation between 
the heavy hole and light hole states, the t-B dependence 
becomes stronger again. 



B. QD's in (111) quantum well 

We now turn to investigate QD's under different strains 
in (111) quantum well. As pointed out before, differing 
from the case of (001) quantum wells, there always ex- 
ists spin mixing between spin- up and -down heavy (light) 
holes in the eigen functions of the Luttinger Hamiltonian 
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FIG. 7: SRT vs. strain at a = 5 nm, rf = 20 nm, S = 1 T and 
T = 4 K. Solid curve: total SRT; Chain curve: SRT induced 
by the hole-phonon scattering due to piezoelectric coupling; 
Dashed curve: SRT induced by the hole-phonon scattering 
due to deformation potential. 



of (lll)-oriented quantum wells. Therefore, hole-phonon 
scattering due to both deformation potential and piezo- 
electric coupling contributes to the spin relaxation even 
there is no strain on QD. Moreover, the effect of strain on 
(lll)-orientcd crystal is also different from the (001) case: 
The strain hardly changes the relative position of energy 
levels of heavy and light holes, but introduces additional 
spin mixing which leads to additional spin relaxation. 

In Fig. 7 we plot the SRT versus in a QD of a = 5 
nm, d = 20 nm and B = 1 T at T = 4 K. It is seen from 
the figure that opposite to the (001) case as shown in 
Figs. 3 and 4, here the SRT decreases rapidly with the in- 
crease of strain in both positive and negative directions. 
This is because that the additional spin mixing intro- 
duced by the strain is the main effect in the present case, 
which makes the scattering rate increases with strain. 
Moreover, one finds that for small strain, the SRT is de- 
termined by the hole-phonon scattering due to deforma- 
tion potential but after > 0.1, hole-phonon scat- 
tering due to the piezoelectric coupling starts to con- 
tribute to the spin relaxation and after |e||/e||| > 0.3 it 
takes over the scattering due to the deformation cou- 
pling and becomes the leading contribution. However 
both contributions should be included in the calculation 
when strain is presented in (111) quantum wells. 

We discuss the diameter and magnetic field depen- 
dence of SRT under different strains. In Fig. 8 the SRT 
is plotted against the QD diameter d for different strains. 
As the SRT with positive and negative strains are sym- 
metrical, we only show the case with e||/e|| > 0. It is 
noted that when e||/e[| = 0.12 the hole-phonon scatter- 
ing due to the deformation potential is dominant but 
when it is 0.28 the scattering due to the piezoelectric 
coupling becomes more important. Similar to the cases 
in the previous sections, the SRT increases with diame- 
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FIG. 8: SRT vs. the diameter d for a QD in (111) quantum 
well at a = 5 nm and B — 1 T under different strains. Curve 
with •: e||/ej} = 0; ■: e||/ejj = 0.12; A: ey/ej} = 0.28. (a): 
T = 0.2 K, (b) T = 4 K. 

ter monotonously when T ~ 0.2 K [Fig. 8(a)] and has a 
maximum when T = 4 K [Fig. 8(b)]. 
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This can be seen in following: In the case when there is 
no inter subband spin mixing and the heavy- and light- 
hole states are separated from each other, then almost all 
the spin mixing comes from the off-diagonal terms of the 
strain Hamiltonian [Eqs. (10) and (11)]. Using the per- 
turbation method, and adopting the lowest four states of 
Hq as basis, the wave functions are written into: 

(r|*T) = (r|0,0,l,^)+C(r|0,0,l,^)+D(r|0,0,l,-^) , 
(r|*^) = (r|0,0,l,-^)+£:(r|0,0,l,^) 

+.F(r|0,0,l,-i) , (34) 

in which C = TnQH/{heBK+Ae), V = niol / {2heBK+Ae), 
£ = -moI*/{2heBK + Ae) and = moH*/{heBK + Ae) 
with Ae = 272?iO — 472^^. H and I are the matrix el- 
ements of strain Hamiltonian [Eqs. (10) and (11)], which 
are independent of the magnetic field. Consequently the 
spin mixing decreases with B. 
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FIG. 10: SRT versus the well width o at d = 20 nm, B = 1 T 
and T = 4 K. Solid curve: (001) quantum well; Dashed curve: 
(111) quantum well. 



FIG. 9: SRT vs. the magnetic field S at a = 5 nm, d = 20 T 
and T = 4 K under different strains. Curve with •: = 0; 

x: e||/4 = 0.01; +: e||/ef| = 0.02; ■: e^^/e^ = 0.12; A: 

In Fig. 9 we show the magnetic field dependence of the 
SRT under different strains. When there is no strain or 
very small strain, similar to the case of QD in (001) quan- 
tum well and our previous investigation on electron spin 
in QD's,^° the SRT decreases with the magnetic field. 
However for a little bigger strain the SRT increases with 
the magnetic field. This is understood that the spin mix- 
ing induced by strain in (111) quantum well makes a 
major contribution to the spin relaxation for sufficient 
big strain. Nevertheless, this mixing decreases with B. 



C. Well width dependence of the SRT 

Finally we investigate the quantum well width depen- 
dence of the SRT of QD's with d = 20 nm and B = 1 T 
at T = 4 K. The SRT's of QD's in (001) and (111) quan- 
tum wells are plotted as function of quantum well width 
a. It is seen that for both cases the SRT's decrease with 
the well width, which is totally opposite to the cases of 
electron spin in QD's^° and quantum wells. This dif- 
ference originates from the fact that for electron spin the 
spin-orbit coupling decreases dramatically with the well 
width. Nevertheless, for hole spin although the spin- 
orbit coupling also decreases with a [see Eqs. (5) and 
(16)], the decrease of the intervals between different en- 
ergy states is faster [see Eqs. (4) and (15)]. Consequently 
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more states are included in the spin-flip scattering chan- 
nel and r decreases with the well width. 



IV. CONCLUSIONS 

In conclusion, we have performed a comprehensive in- 
vestigation on hole spin relaxation in GaAs QD's con- 
fined in quantum wells along (001) and (111) directions 
by exactly diagonalizing the hole Luttinger Hamiltonian. 

We find for QD's in (001) quantum wells with small 
well width where only the lowest subband is involved, the 
SRT increases with the QD diameter at very low temper- 
ature {e.g., 0.2 K) or first increases until it reaches a max- 
imum and then decreases at higher temperature (e.g., 4 
K). These features are opposite to those of electron spin 
in QD's. Moreover, unlike the case of electron spin where 
the SRT is mainly determined by the electron-phonon 
scattering due to the piezoelectric coupling, here only 
the hole-phonon scattering due to the deformation po- 
tential contributes to the spin relaxation. Strain changes 
the relative positions of energy levels of heavy hole and 
light hole. A positive strain increases the energy gap be- 
tween heavy hole and light hole and enhances the SRT. 
A negative strain decreases the gap and reduces the SRT 
until the interchange of the lowest energy states from 
heavy hole to light hole. After that the SRT increases 
again. Moreover, very close to the interchange point, as 
the energy levels of heavy hole and light hole are very 
close to each other, the hole-phonon scattering due to 
piezoelectric coupling contributes to the spin relaxation 
too. For large well width where multi-subband effect is 
important, hole-phonon scattering due to the piezoelec- 



tric coupling contributes to the spin-flip scattering with 
or without strain. Nevertheless the SRT is still mainly 
determined by the scattering due to deformation poten- 
tial except at the interchange point. The magnetic field 
dependence of the SRT is also discussed. 

For QD's in (111) quantum well things are quite dif- 
ferent from those for QD's in (001) quantum well: Hole- 
phonon scattering due to both piezoelectric coupling and 
deformation potential contributes to the spin relaxation 
and should be both included in the calculation, regardless 
of the well width; Strains can hardly change the relative 
positions of energy levels of heavy hole and light hole 
but introduce additional spin mixing. Therefore the SRT 
decreases rapidly with strain; The SRT decreases with 
magnetic field like the case of QD in (001) quantum well 
and our previous investigation on electron spin in QD's 
when there is no strain or very small strain. However for 
strain which is big enough that the spin mixing is mainly 
determined by it, the SRT increases with B. 

Finally we show that the hole SRT decreases with well 
width for QD's in both (001) and (111) quantum wells 
which is totally opposite to the cases of electron spin in 
QD's and quantum wells. 
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